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^> \ Abstract IRAS22475+5939 has been well researched by previous astronomers. But 



we still get some new characteristics about it, using the first observations in lines of 
CO J = 2 - l, 13 COJ = 2 - l, 13 COJ = 3 - 2 by the KOSMA 3 m telescope. 
The mapping of the intensity ratio of,Eiroa 1981,Felli81,Eiroa 1981 13 COJ = 3-2 
and 13 CO J = 2 — 1 shows the distribution of the temperature with two peaks, which 



don't coincide with IRAS22475+5939 source and the center of the HII region, but at 
the edge of the HII region. The overlays of the Spitzer IRAC 8 pirn and CO contours 
indicate that they are associated with each other and the strongest polycyclic aromatic 
hydrocarbons (PAHs) emission is at the position of IRAS22475+5939 source. While the 
IRAS LRS spectrum at 7 pm ~ 23 pm and the PHT-s spectrum at 2 pm ~ 12 pm of 
IRAS22475+5939 source also exhibit strong PAHs emission characters at the main PAH 
bands. The diversity of PAH family should be responsible for the plateaus of PAHs emis- 
sion in the PHT-s spectrum and the IRAS-LRS spectrum. An analysis and modeling in 
infrared bands suggest that IRAS22475+5939 is more likely to be a Class I YSO Where 
ly-v this is the case, the star is likely to have a temperature Teff ~ 9995.8 K, mass ~ 

15.34 M Q ,lumilosity ~ 1.54 x 10 4 L Q , andage ~ 1.54 x 10 4 yr. The model shows 
that the circumstellar disc emission is important for the wavelength between 1 and 10 pm, 
otherwise, envelope fluxes for A > 10pm. The bipolar outflow is confirmed in the molec- 
ular cloud. The excited star of the HII region has the chance to be the driving source of 
the outflow. The high resolution is required. 
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1 INTRODUCTION 

Researches of massive star formation have been getting significant attention in the recent years. But 
the mechanism of how massive stars form is still debatable, whether it is through the process similar to 
that for the low mass stars, ie., through disk accretion and driving molecular outflows (Shu et al. 1987), 
or through alternative scenarios-coalescence of low mass stars (Bonnell et al. 1998). But the recent 
observations show that outflows are common in the massive star formation regions (Zhang et al. 2001) 
and it is almost clear that stars at least up to late-O spectral types form primarily through disk accretion 
(Varricatt 2012). HII regions will form when the massive stars ionize their surroundings and provide 
information about massive star formation within molecular clouds (Heyer et al. 1989;Churchwell 2002). 
Therefore, the forming stars will have a big influence on their surrounding environment during the 
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evolutional processes. It is an indirect way to study the properties in the star formation regions in order 
to know more about how massive stars form. 

IRAS22475+5939 has been detected in the sharpness HII region S 146 (Zuckerman & Evans 1974). 
The distance of IRAS22475+5939 from us is about 4.7 kpc (Henkel 1986). Eiroa (1981) believed S146 
was excited by an IRS 1, locating at position a(B1950) = 22 h 47 m 29 s .7, (5(51950) = +59°38'55" ± 7 
at 2.2 nm, a(£1950) = 22 h 47 m 29 s .5, 5(51950) = +59°39'01" ±3 at 0.9 nm. Blair (1978) 
found a H 2 maser associated with IRAS22475+5939. A 2MASS source 22492900+59545600( a = 
22 h 49 m 29 s .06, 8 = +59°54'55".7, J2000 ) was detected by Wang(1997). Felli (1981) showed a bipolar 
outflow with the optical and 6 cm emissions. The bipolar molecular outflow was also confirmed from 
surveys of the CO J = 1 - (Yang & Wu 1998), CO J = 2 - 1 (Jiang et al. 2001 ) and CO J = 3 - 2 
(Wu et al. 2005) and Jiang (2001) thought a massive star was forming in the molecular cloud. Guan 
(2008) made a mapping survey of the massive CO cores with the CO, 13 CO, 18 CO J = 1—0 lines. 
The contour of 13 CO J = 1 — was studied to be associated with MSX 8 /im emission (Guan & Wu 
2008). The IRAS-LRS (Low Resolution Spectra) spectrum was analyzed by Muizon (1990), Volk (1991) 
and Chen (1995), respectively. The LRS spectrum shows strong [NeII](12.8/zm), [NeHI](15.5/z?7i), 
[SIII](18.7/xm), PAHs emission (7.7 /zm, 8.6 fxm, 11.3 /im) as well as the sillicate absorption feature at 
9.7 [im (Muizon et al. 1990; Volk & Kwok 1991; Chen et al. 1995). 

In this paper, we show the existence of the outflow and discuss the properties in the observed region. 
13 CO is optically thin and can trace the internal region of the molecular core. And the intensity ratio 
Ri 32/21 (ratio of 13 CO J = 3 — 2 and 13 CO J = 2 — 1) contains the information about the temperature 
distribution in the molecular cloud, and mid-infrared emission at 8 um is thought to be from small dust 
grains and polycyclic aromatic hydrocarbons (PAHs), and excited by the UV radiation leaking from the 
HII regions (Leger & Puget 1984; Deharveng et al. 2003,2005). We also use the Stokes I image from the 
observations for the 1.4 GHZ NRAO VLA Sky Survey(NVSS) to trace the HII region in order to know 
the relationship between the molecular cloud, the intensity ratio and the HII region. We describe the 
observation in section 2. The results and discussion of the properties are in section 3. section 4 makes a 
conclusion. 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 The observational data 

We carried out observations toward IRAS22475+5939 ( a = 22 h 49 m 29 s .4, 8 = +59°54'54".00, 
J2000) source in CO J = 2 - l, 13 CO J = 2 - l, 13 CO J = 3 - 2 lines using the KOSMA 3 m 
telescope at Gomergrat, Switzerland in April 2004. The half-power beam widths of the telescope at ob- 
serving frequencies 230.538 GHZ, 220.399 GHZ, 330.588 GHZ are 130", 130" and 80" , respectively. 
The pointing and tracking accuracy is better than 10". The DSB receiver noise temperature is about 120 
K. The medium and variable resolution acousto-optical spectrometers have 1501 and 1601 channels. 
The channel widths of 248 MHZ and 544 MHZ correspond to velocity resolutions of 0.21 km • s _1 and 
0.29 km • s"\ respectively. The beam efficiency is 0.68 at 230 GHZ and 220 GHZ, but 0.72 at 330 GHZ 
and 345 GHZ. The forward efficiency is 0.93. The 80" resolution of the J = 3 — 2 data was convolved 
to 130" with an effective beam size -\/(130 2 — 80 2 ) = 102". The correction for the line intensities 
to the main beam temperature scale was made using the formula T m b = (F c g/B g * T^). The data 
were reduced by the software CLASS (Continuum and Line Analysis Single-Disk Software) and GREG 
(Grenoble Graphic). 

2.2 Archival data 

The Spitzer IRAC 8 um imaging is available toward the direction of IRAS22475+5939 at an angular 
resolution of ~ 2". S146 was detected well by the 1.4 GHZ NRAO VLA Sky Survey(NVSS). Its 
resolution is 45" and a limiting peak source brightness is about 2.5 mJy/beam. We have compiled the 
SED for IRAS22475+5939 source using near-IR (JHK) flux from the 2MASS All-Sky Point Source 
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Fig. 1 Left: spectra of CO (J = 2 - 1), 13 CO (J = 2 - 1) and 13 CO (J = 3 - 2) at the 
position of (-1, 0). Right: channel maps of the corresponding isotopes 



Catalog (PSC), mid-IR MSX four bands from the Midcourse Space Experiment (MSXC6), and the 
IRAS data from IRAS Point Source Catalog v2.1 (PSC). The IRAS LRS spectrum is obtained from the 
LRS database at the University of CalgarjQ (Hodge et al. 2004;Kwoket al.1997). The LRS was a slitless 
spectrometer and well-suited to point sources whereas confusion arises in the case of extended source. 
Two wavelength channels are recorded simultaneously: Band 1 (7.7 — 13.4 fim) and Band 2 (11.0 — 
22. 6 lira). The wavelength resolution A/AA varies between 10 and 60; it increases with wavelength 
inside each wavelength band, the resolution being systematically somewhat lower in band 2 than in band 
1. PHT-s consist of a dual grating spectrometer with two 64-element arrays that span the 2.5 — 4.9 jim 



See http://www.astro.wisc.edu/protostars 
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and 5.8—11.6 \im spectra regions. The Auto Analysis Result produced by the Off-Line Processing(OLP) 
from the ISO Data Archive(IDA) is available (Hodge et al. 2004;Lemke et al. 1996) 

3 RESULTS AND DISCUSSION 
3.1 The outflow 

3.1.1 The molecular spectral lines 

Figure 1 shows the spectra and channel maps of CO J = 2 - 1, 13 CO J = 2 - 1, 13 CO J = 3 - 2. The 
spectra are at the positions (-1,0). Each spectrum shows a broad width as well as the almost same Vslr 
and shape. The showing spectra are the strongest for 13 CO J = 2 — 1, 13 CO J = 3 — 2. The optically 
thick line CO J = 2 — 1 is not gaussian shape and has wings. The blue wing and the red wing are not 
asymmetrical. Table 1 presents the observational parameters of IRAS22475+5939 source. The ranges 
of the FW come from P- V diagram, which shows the high-velocity gas exists. 

Table 1 Observation physical parameters of IRAS22475+5939 
Source 

Name T mb FWHM FW Vlsr 
(K) (km-s- 1 ) (km-s' 1 ) (km ■ s" 1 ) 

CO(J = 2-l) 11.20 4.99 ±0.03 12.37 -49.60 ± 0.01 

13 CO(J = 2-l) 3.02 3.90 ±0.08 6.33 -49.80 ± 0.04 

13 CO(J = 3-2) 2.55 4.79 ±0.17 7.87 -49.69 ± 0.07 



3.1.2 The outflow and its physical parameters 

Table 2 Physical parameters of the outflow 



Name 
IRAS 


Wing- 


N(H 2 ) 
(xl0 20 cm- 2 ) 


M 
(M ) 


t d 
(xl0 5 yr) 


M 

(xlO^Moyr- 1 ) 


22475 ± 5939 


blue 
red 


1.00 
0.15 


278.79 
26.63 


5.37 
3.25 


3.86 
0.80 


Name 
IRAS 


Wing 


F 
(M Q kms _1 • yr" 


P 

- 1 ) (Mekms- 1 ) 


E LmbcIi 

(xl0 46 erg) (L ) 


22475 + 5939 


blue 
red 


1.93 x 10 -3 
3.99 x 10" 4 




1028.74 
129.699 


7.59 1.17 
1.26 1.23 



We show the outflow and its corresponding P-V diagram in Figure 2. The integrated ranges of the 
blue wing and the red wing are determined by the P-V diagram and gaussian fit. We find they overlap 
with each other, which is likely they are in the direction of the sight or the resolution of the telescope 
is too low to identify the direction of the outflows. And IRAS22475+5939 as well as its associated 
H2O maser are almost in the center of the wings. Maybe IRAS22475+5939 is the driving source of 
the outflow. The physical parameters are presented in Table 2. Assuming it is local thermodynamic 
equilibrium (short for LTE latter) and CO in the outflow is optically thin, we can calculate the col- 
umn density of the outflow using the below formula (Scoville et al. 1986) under the assumption of 
N(CO)/N(H 2 ) « 1CT 4 (Dickmanetal.1978): 
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Fig. 2 Left panel: outflow contours of CO (J = 2 — 1), the integrated ranges in blue wing 



and red wing are —55.78 km • s 



-50.80 km • s~\ -46.05 km • s" 



-43.41km- 



s , respectively. The contour levers are 30% ~ 90% of each wing's peak value. The star 
is IRAS22475+5939, the triangle is the H 2 maser, and the square is the 2MASS source 
22492900+59545600. Right panel: P-V diagram. The contour levers are 0.7 - 3 by 0.6 K, 
3 ~ 12 by 1.0 K, 12 ~ 20 by 1.5 K. The vertical lines indicate the beginning of the blue and 
red wings, respectively. 



N = 10 5 x 



3k 2 



4h7T 3 ^ 2 I/' 



hv J T cx 

exp( ) 



mv/6k(J + 1) 



2kT cx exp(-h^/kT ox ) 



1-e" 



T mb dv 



(1) 



where T ox is the excited temperature, which is got from the equ(9). J is the lower lever of the transition. 
/i = 0.112D. 

The other physical parameters coming from the below formulas (Xu & Wang 2010) : 



M = um H2 SN(H 2 )/2x 10 

P = MV 

E = MV 2 
t d = R/V 

F = P/t d 
M = P/t d v w 

LMcch = E/td 



:>,:! 



(2) 

(3) 

(4) 
(5) 
(6) 

(7) 
(8) 



where M is the outflow's mass, td, P, E, V, R, F, LMech are the dynamic time, the momentum, the 
energy, the mean velocity of the gas relative to the Vlsr> the size of the wings, the driving force, the 
mechanical luminosity. The mean atomic weight of the gas u = 1.36, S is the area of the outflow. 
From these parameters, we can see the outflow is massive and energetic than the low-mass stars, so 
we further validate that it is possible to have massive stars forming in this region. At the same time, 
through the outflow and its parameters, we find the blue outflow is far larger than the red one, and both 
are elongated from east to west, very similar to the 13 CO J = 3 — 2 molecular core that is presented 
latter. The explanation is that the surrounding gas might mix with the blue wing. But we can not make 
sure because of the low resolution of the telescope. 
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Fig. 3 Integrated intensity diagrams of the core emission. In each map, the integrated range 



is from -50.80 km • s _1 to -46.05 km ■ 
value. 



s : . The contour levels are 30% to 90% of the peak 



3.2 The molecular core 



Figure 3 shows the structure of the molecular core. The integrated range is determined by the P-V 
diagram and the gaussian fit due to the relatively large noise. It is more accurate than one determined 
by each of them. In the diagrams, we find the core structure from 13 CO J = 3 — 2 is elongated from 
east to west, but almost round from the transition CO J = 2 - 1, 13 CO ,1 = 2-1. IRAS 22475+5939 
and H2 O maser coincide with the peak of cores. We also calculate the core's physical parameters by the 
expressions (Garden et al. 1991) under the assumption of LTE: 



T« ( co) = £(m(i + £[^ 



hi//k 



exp(hi//kTb g ) — 1 



L )y 



(9) 



r( 13 CO) = -ln{l- 



kT 



nib 1 



hv cxp(h^/kT C x) — 1 cxp(hi//kTb g ) — 1 



} 



(10) 



where Tb g = 2.732 K is the temperature of the cosmic background radiation, f is the beam filling factor, 
here we assuming / = 1. Assuming t(CO)/t( 13 CO) = [CO]/[ 13 CO] = 89 (Langct al.1980). We can 
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use another method to calculate the optical depths: 

T mb (CO) 1 



exp[-r(CO)] 



(ID 



T mb pCO) l-cxp[-r(i3CO)] 

The two calculated values are summarized in Table 3. Comparing with them, we find they are the same 
in the error scale. Therefore we can draw a conclusion that our assumptions are right and the abundance 
ratio [CO]/[ 13 CO] in this region is almost the same as in our solar system. The core's physical parame- 
ters are in Table 4. From the figures and tables, we can know that the core is isolated and massive, which 
can provide the environment for the high-mass star formation. 

Table 3 Calculated results of the optical depth by two 
methods 



Name 



T\ 



Tl 



i CO(J 



0.28 



0.28 



Notes: n is derived from the fist method, T2 is from the second method. 



Table 4 Physical parameters of the core 



Name 


J- ex 

(K) 


T 


/y(COJ = 2-1) 
(xl0 18 cm- 2 ) 


N(R 2 ) 
(xl0 22 cm- 2 ) 


M 
(xlO 3 M ) 


IRAS22475 + 5939 


17.35 


24.92 


1.03 


1.03 


3.56 



3.3 The intensity ratio Ri 
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Fig. 4 Grayscale shows the 1.4 GHZ NVSS image. The intensity ratio Ri 32/21 is the green 
contour, the blue and red contours are the integrate intensity maps of 13 CO (J = 2 — 1), 
13 CO (J = 3 - 2), respectively. The star is IRAS22475+5939, the triangle is the H 2 maser, 
and the square is the 2MASS source 22492900+59545600. 



Figure 4 gives the contours of 3 CO and the intensity ratio Ri 32/21 overlapping on the 1.4 GHZ 
NVSS grayscale diagram which can trace the structure of the radio source, here we can identify the 
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Fig. 5 The molecular cores from the transitions CO(J = 2 — 1), 13 CO (J = 2 — 1), 
13 CO (J = 3 - 2) (top left, top right, bottom left) and the outflow of CO (J = 2 - 1) 
emission (bottom right) overlaid on the grayscale Spitzer 8 um emission. The synthesized 
beam size (~ 1') is shown in the lower left-hand corner of the each panel of KOSMA 3m 
telescope. 



radio source is a HII region (Lockman et al. 1989). The intensity ratio indicates the gas temperature 
distribution in the region (Qin et al. 2008). From Fig 4, we can find the morphology of the intensity ratio 
map is similar to a triangle and has three peaks. The peaks are not in the center of the molecular cloud 
core, the biggest peak is in the northwest, about 2' away from the center. The value of it is ~2.12 > 1, 
which is bigger than the value (< 1) detected in Cepheus B (Beuther et al. 2000). Maybe it is related with 
HII regions (Wilson et al. 1997). The second peak and the third peak are in the direction of southwest 
and east, respectively. All of three are at the edge of the HII region. This shows the temperature at the 
edge is higher than the center. Maybe the hot gas of the molecular cloud was pushed away to the edge 
and heated by the HII region. This explanation can be supported by the distribution of the HII region and 
the molecular core. Additionally, we find 13 CO higher J can trace the warmer region. And the center of 
the molecular core coincides with that of the HII region. 
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Fig. 6 Left panel: the IRAS-LRS spectrum coming from LRS database at the university of 
Calgary, consistent of two wavelength band: 7. 7— 13.4 fim, 11.0 — 22.6 \im. Right panel: the 
PHT-s spectrum lying in the 2 — 12 /j,m range. 

3.4 polycyclic aromatic hydrocarbons (PAHs) 

Fig 5 shows the superposition diagrams of the CO isotopes' contours and CO outflow on the Spitzer 
8 fim emission. The Spitzer IRAC 8 fim emission is primarily due to the 7.7 urn and 8.6 fim PAHs 
features. The grayscale diagram shows the distribution of PAHs in the HII region. We can see that the 
strongest PAHs emission coincides with IRAS22475+5939 source, which shows the PAHs emission 
probably is excited by the UV radiation from IRAS22475+5939. 

Fig 6 shows the LRS spectrum and the PHT-s spectrum, spanning the wavelength 7.7 — 23 fim and 
2 — 12 fim, respectively. They are dominated by PAHs emission features at the so-called unidentified in- 
frared bands (3.3,6.2, 7.7,8.6, 11.2 fim) and weaker bands0(Peeters 2011). The PHT-s spectrum shows 
the narrow and strongest emission at 7.6 fim at the wavelength range 2 — 12 fim and an absorption at 
9.7 fim, but the strongest emission for the LRS is at 12.0 fim at band 1 . It is more believable for the PHT- 
s due to its higher resolution than the results of the LRS. And a broad emission plateau at 16.4—17.4 fim 
maybe resolute from the low resolution or the diverse PAH family. Muizon (1990) pointed out that the 
prominent emissions at 12.8, 15.55, 18.7 fim are [Nell], [Nelll], [SIII] lines,respectively. 

3.5 The spectral energy distribution 

Table 5 Results from SED fit- 
ting 



Parameter 


SM 
(M ) 


age 

(yr) 


Teff 
(K) 


DM 

(M ) 


Disk accretion rate 

(Moyr" 1 ) 


Luminosity 

L 


A v 


Best fit values 


15.34 


1.54 x 10 4 


9995.8 


0.13 


1.56 x 10"° 


1.54 x 10 4 


4.53 



Notes: SM is short for stellar mass, DM is for disk mass, the luminosity is the total luminosity of the star, 
A v is the foreground extinction. The best fit values are derived the average values for the best top 10 models. 

The 2MASS, MSX, and IRAS data are used to construct the SED of the IRAS22475+5939 source. 
Color correction is applied to the IRAS data using the correction factors given in the point source 



/ This includes bands at 3.4, 3.5, 5.25, 5.75, 6.0, 6.6, 6.9, 7.2 - 7.4, 8.2, 10.5, 10.8, 11.0, 12.0, 13.5, 14.2, 15.8, 16.4, 
16.6, 17.0, 17.4, 17.8, 19.0 /im 
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Fig. 7 Left paneLthe SED of the IRAS 22475+5939 source. The filled circles show the data 
from 2MASS (JHK), MSX (ACDE), IRAS (12, 25, 60, 100 urn). The continuous line shows 
the best fit model and the grey lines show subsequent good fits for (% 2 — xlest fit ) P er c ' ata 
point < 3. The dashed line corresponds to the stellar photosphere for the central source of the 
best fitting model, as it would look in the absence of the circumstellar (but including interstel- 
lar extinction). Right panel: various emission components making up this model, details in the 
text. Suffice to remark that the circumstellar disc (green curve) is important for wavelengths 
between 1 and 10 /im. 



catalogue (Beichman et al. 1988). The SED fitting tool of Robitaille et al. (2007) is available on-line to 
model the SED. The SED plot is showed in Fig 7 and its fitting parameters in Table 5. We derived a mass, 
luminosity and temperature of 15.34 M , 1.54 x 10 4 L Q , 9995.8 K respectively from the fit. And the 
average foreground extinction is ~ 4.53. From the SED and the active accreting mass, IRAS22475+5939 
is probably to be a class I protostar. 

We have also indicated the differing component of flux which make up the model YSO spectrum 
(right panel of Fig 7). In this case, the total flux is indicated as black, the stellar flux as blue, the stellar 
photospheric flux is showed in the dashed line (this is the flux prior to reddening by circumstellar dust), 
the disc flux as green, the scattered flux as yellow, the envelope flux as red and the thermal flux as 
orange. Unless otherwise stated, the results include the effects of circumstellar extinction, but not of IS 
extinction. The also assume a representative distance of 1 kpc. 

It is apparent, from the latter modeling, that the central star and disc components of emission are 
responsible for the MIR/NIR emission, but envelope fluxes are of most importance at A > 10 [im. 



3.6 Discussion 

The morphology of the cores are firstly showed from the lines CO J = 2 — 1, 13 COJ = 
2 — 1, 13 COJ = 3 — 2, as well as the corresponding parameters. Comparing with the re- 
sults in CO, 13 CO, 18 CO J = 1 - (Guan et al. 2008), the V LS R and the derived physic pa- 
rameters are litter difference. And the structure of 13 COJ = 1 — detected by them is simi- 
lar to 13 CO J = 2 — 1. The outflow traced by CO J = 2 — 1 is more extended than that of 
CO J = 3-2 (Jiang et al. 2001). Yang(1998) only analyzed the IRAS data and obtained that 
IRAS22475+5939 may be a massive YSO with high luminosity. Due to the SED in infrared bands, we 
suggest that IRAS22475+5939 is likely to be a class I YSO with mass 15.34 M Q ,lumnosity 1.54 x 
10 4 L Q , age 1.54 x 10 4 vt,Teff 9995.8 K. Its foreground extinction is about 4.53 mag. Muizon 
(1990) suggests the PAH characteristics at 7.7, 8.6, 11.3 /im, silicate absorption features at 9.7 /j,m 
and the prominent emissions at 12.8/jm([NcII]), 15.55 /im([NeIII]), 18.7/im([SIII]) according to 
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the IRAS-LRS spectrum. But we combine the IRAS-LRS spectrum with the PHT-s spectrum, find- 
ing the PAHs emissions are showed at all the PAH bands, especially at the main PAH bands 
(3.3, 6.2, 7.7, 8.6, 11.2 /im) and 16.4 — 17.4 [im. This indicates the diverse PAH family in the molecular 
cloud. The HII region is likely to be excited by IRl-an 07V-07.5V star (Eiroa et al.1981). Its position 
almost coincides with IRAS22475+5939. Therefore, the star has the chance to drive the outflow. It is 
necessary to have high resolution observation to identify their relationship. 

4 CONCLUSION 

We observed the IRAS22475+5939 source using three spectral lines CO J = 2 - 1, 13 CO J = 2 - 1, 
13 CO J = 3 — 2. The Vslr of them are almost the same. The massive and energetic outflow is further 
verified by us, showing that this region is a high-mass star formation region. The integrated intensity 
maps of the cores tell us the molecular cloud is isolated. And the core's mass is larger than the low-mass 
cores. The intensity ratio Ri 3 , /21 map indicates the gas temperature varies at different positions and the 
maximum value is bigger than 1, which is larger than the regions without HII regions (Wilson et al. 
1997). And the peaks of the intensity ratio are at the edge of the HII region. Maybe the HII region blows 
away the hot molecular gas from the center to the edge and heats the gas. This causes the temperature is 
higher at the edge than the center of the molecular core. And 13 CO higher J transition traces the warm 
regions. PAH features are demonstrated in the Spizter IRAS 8 [im emission and the spectra of IRAS- 
LRS and PHT-s, suggesting a rich PAH molecules in the cloud. IRAS22475+5939 may be a massive 
and luminous class I protostar with a active disc. But considering the position of the excited star of the 
HII region, the excited star is able to drive the outflow. Maybe it is the driving source of the outflow. The 
high resolution observations is needed for better results. 
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